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Abstract 

(N ■ 
£> J We develope the method of improvement of perturbative theory in QCD, applied 

to any polarization operator. The case of polarization operator Ii(q 2 ), corresponding 
to the process e + e~ — > hadrons is considered in details. Using the analytical prop- 
erties of II((/ 2 ) and perturbative expansion of H{q 2 ) at q 2 < 0, ImU(q 2 ) at q 2 > 
is determined in such a way, that the infared pole is eliminated. The convergence of 
q^ . perturbative series for R(q 2 ) = a(e + e~ — ► hadrons) /(e + e~ — > fi + ^~) is improved. 

After substitution of R(q 2 ) into dispersion relation the improved Adler function D(q 2 ) 
with no infrared pole and frozen a s (g 2 ) has been obtained. A good agreement with 
J* ■ experiment has been achieved. 

<D : 

43 ■ PACS: 12.38. Cy 

>: 

It is well known for years, that perturbative calculations of amplitudes in field theories are 
legitimate, if virtualities of all external legs k 2 are negative, k 2 < . In this case polarization 
operators and vertex functions (two- and three-point functions) are off mass shell and have 
no singularities. However, in order to get physical predictions, it is necessary to go to k 2 > 
(at least for some of them). This can be achieved by analytical continuation using known 
analytical properties of amplitudes. A typical example is e + e~-annihilation into hadrons in 

QCD. The total cross section a(e + e > hadrons) is proportional to imaginary part of the 

polarization operator li^q 2 ) of virtual photon at q 2 > 0. At large enough q 2 < H^q 2 ) 
can be calculated perturbatively in QCD in terms of the expansion over the running coupling 
constant a s (q 2 ) = Air/ (3oln(— q 2 /A 2 ). The tensor structure of II^(g) is 

IV(g) = (g M g„ - g 2 ^)n(g 2 ) (1) 

and II(g 2 ) is an analytical function of q 2 in the whole q 2 complex plane with a cut along 
positive q 2 semiaxes. Analytical continuation from q 2 < to q 2 > results in substitution 
of ln(q 2 /A 2 ) — ► ln(Q 2 /A 2 ) —in, Q 2 —\ q 2 \. Since small a s corresponds to large ln(Q 2 /A 2 ), 
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the standard procedure (see, e.g. |I[) is to consider ln(Q 2 /A 2 ) as large comparing with 
ir and to perform the expansion in 5 = tt / ln(Q 2 / A 2 ) . However, practically it is not a 
good expansion parameter. At typical scale of e + e~-annihilation, Q 2 ~ 10 GeV 2 and A ~ 
300 - 400 MeV, ln(Q 2 /A 2 ) « 4 - 5 and S « 0.7. 

In this paper we present the systematic method of improvement of perturbation theory 
in QCD, which is free from this drawback. Besides e + e~ annihilation, this method can 
be applied to any polarization operators, for example, to those used in the QCD sum rule 
approach. The idea of the method has been suggested by Radyushkin 0] and considered also 
by Pivovarov ||, but many important features and consequences of this method were not 
touched in f2|,[|| (particularely, eq. (|To| ) below was not obtained and analysed). Consider 
the Adler function 
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where R(s) 
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hadrons)/a(e + e — > fi + fi ). In the parton model R{s) = R p 



3 J2 e 2 where e q is the charge of the quark of flavour q. It is convenient to write: 

D(q 2 ) = R p (l + d(q 2 )), R(q 2 ) = R p {l + r(q 2 )), U(q 2 ) = R p (l + p(q 2 )) 
From (Q), there follows the equation: 

d(q ) = -q 



dq 1 



The solution of (|4T) is 
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r(q 2 ) is proportional to discontinuity of p(q 2 ) at q 2 > 0: 
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At negative q < 0, q — —Q the perturbative expansion of d(Q ) in MS renormalization 
scheme is known up to the third order 0,0: 

d(Q 2 ) = (1 + d x a + d 2 a 2 ), d x = 1.986 - 0.115N f 

d 2 = 18.244 - 4.216A^/ + 0.086A^J (7) 

where a(Q 2 ) = a s (Q 2 )/n, Nf is the number of flavours and the small gluon-gluon scattering 
terms are omitted. With the same accuracy, three loops expression for a s (Q 2 ) in MS is 
given by @, 0| 
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Substitution of (0), (§) into (|5|) and (||) leads to perturbative corrections up to the third order 
in the physically measurable quantity r(g 2 ).By taking the discontinuity, any dependence on 
the normalization point a 2 is eliminated. It should be stressed, that in such calculation 
no expansion in n / ln(Q 2 / A 2 ) is performed, the only assumptions used(in fact, they are not 
assumptions, but theorems), are the analytical properties of II(g 2 ) . The result is (q 2 > 0): 
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Essential features of (10 ) are: 1) unlike standard perturbation expansion r(q 2 ) has no infra- 
red poles at q 2 = A 2 and tends to finite limit at q 2 — > 0, r(0) = 4//3 ~ 0.414, corresponding 
to a s (0) = 47t//3o ~ 1-3, all higher order terms vanish; 2) the convergence of perturbation 
series in (|10"D is much better, than of the standard ones: the ratio of the second order term 
to the first order one is smaller than 0.2 everywhere and the ratio of the third order term 
to the second one is smaller than 0.1, while in the standard approach the ratio of the third 
order term to the second order one exceeds 0.5 below q 2 = 10 GeV 2 and is larger than 1 at 
q 2 < 1 GeV 2 ; 3) the a s corrections, given by (|T0"D, are remarkably smaller in the low energy 
domain, q 2 < 5 GeV 2 , than the standard ones, e.g., at q 2 = lGeV 2 r(eq.l0)/r stand = 0.72. 

In order to get r(q 2 ) at some lower values of q 2 , starting from the value of a s at Z-boson 
mass, a s (m z ), which is rather well known now, it is necessary to use renormalization group 
equations ([8]) and perform matching at the thresholds of new flavours production b and c 
quarks and, if we would like to go to q 2 < 1 GeV 2 , even at s-quark threshold. This matching 
procedure introduces some uncertainty. The matching may be performed at 2m q (or, what is 
practically equivalent, at m-f and mj/f) or at m q . The former seems to be preferable in the 
case of e + e~-annihilation for the evident reasons. There are some arguments in the favour of 
the latter choice, based on minimal sensitivity of the results to the matching point value || . 
In the standard approach the results are rather sensitive to the choice of matching points 
- the ratio of r's in the two above mentioned cases is about 1.10 — 1.15 at q 2 < 5 GeV 2 . 



We calculated the q 2 dependence of r(q 2 ) starting from the point g 2 



m: 



and going down. 



The value r(m z ) was found from the requirement, that a s {m z ) determined in the standard 



way is equal a s (m 2 z ) = 0.119 ± 0.002 @. Then A 5 = 23011? MeV was found 



evolution down to lower energies the matching of r(q ) at the masses of T, 
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J/ijj and if was 



performed resulting in A 4 = 
found particularly, that nr( 



33511 M eV, A 3 = 4141^ MeV, A 2 = 490+1} MeV. It was 
1 GeV 2 ) = 0.41. If, instead of matching at quarkonium masses 



the matching at m q would be done, the values of r(q 2 ) would be only 4% smaller below 
q 2 = 10 GeV 2 practically independent of Q 2 . Therefore, in this aspect this method has 
also some advantages. The performed comparison with experiment demonstrated a good 
agreement, starting from ^/q 1 = 0.7 GeV, i.e. in much broader interval, than in recent paper 

After substituting r(q 2 ), given by ( [TDD ^° dispersion relation (^j) we get back Adler func- 
tion D(q 2 ). Obtained in this way improved D(q 2 ) impr has correct analytical properties and 
no unphysical singularities. Therefore, the adopted here procedure has a serious advantage: 
the required analytical properties of Adler function are restored. If D(q 2 ) impr is represented 
in terms of improved effective QCD coupling constant a s (q 2 )i mpr , this would mean, that 
ot s (q 2 )impr has no infrared pole and is frozen at q 2 — > 0. In this aspect our approach has 
some resemblance to the approach by Shirkov and Solovtsov JTTJ (see also ||12||), where the 



condition of analyticity of a s (q 2 ) in the cutted q 2 complex plane was imposed. The difference 
is, that we exploit the analyticity of polarization operator II(g 2 ), which is a strict result in 
field theory and do not use any hypothesis about analyticity of a s (q 2 ). 

The method to define the improved QCD coupling constant a s (q 2 )i mpr through D(q 2 ) impr 
looks very promisable. For this goal, perhaps, the most suitable is to use not MS, but 
the Brodsky, Lepage, Mackenzie renormalization scheme |TJ|, where e + e~-annihilation is 



considered as a basing process for definition of a s (q 2 ) with no higher order a s corrections. 
This problem requires further investigation. 

The presented above method can be applied to the treatment of perturbative corrections 
to any polarization operators and, therefore may result to improvement of QCD sum rule 
approach. 
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